Hahn-Tsai复合材料的非线性杂交应力有限元方法

成功建立了Hahn-Tsai复合材料模型的非线性杂交应力有限元方程,采用Newton-Raphson迭代法求解结构的非线性位移方程。在迭代过程中,为了提高计算效率可采用简单迭代法由节点位移求解单元应力场。但是,当载荷增加到一定程度以后,非线性应力场由于循环迭代而无法收敛,显然,一般的加速方法不能解决这种循环迭代的发散问题。因此,本文发展了一种确实有效的非线性应力场迭代新方法,在不增加计算工作量的情况下,不仅极大地提高了收敛速度,而且对于较大载荷也能够很好地收敛,从而解决了大载荷下非线性杂交元方法失败的关键问题。数值算例表明该方法是确实可行的。     

Dynamics of milling processes with variable time delays

A milling-process model with a variable time delay associated with each cutting tooth is presented in this article. The source of this variable time delay is the feed rate. The effect of the feed motion on the entry cutting angle, the exit cutting angle, and the amplitude of feed mark is also discussed. Loss-of-contact effects are also considered. The system dynamics is described by a set of delay differential equations with periodic coefficients and variable time delays. A semi-discretization scheme is presented for analyzing the stability of periodic orbits of this system. The analysis provides evidence of period-doubling bifurcations and secondary Hopf bifurcations. Good agreement is found between the numerical results obtained from this work and the results of related experimental studies.     

Measurement of energy release rate and energy flux of rapidly bifurcating crack in Homalite 100 and Araldite B by high-speed holographic microscopy

High-speed holographic microscopy is applied to take three successive photographs of fast propagating cracks in Homalite 100 or in Araldite B at the moment of bifurcation. Crack speed at bifurcation is about 540 m/s on Homalite 100, and about 450 m/s on Araldite B. From the photographs, crack speeds immediately before and after bifurcation are obtained, and it is found that discontinuous change of crack speed does not exist at the moment of bifurcation in the case of Homalite 100, but exists in the case of Araldite B. From the photographs, crack opening displacement (COD) is also measured along the cracks as a function of distance r from the crack tips. The measurement results show that the CODs are proportional to √r before bifurcation. After bifurcation, the CODs of mother cracks are proportional to √r, though the CODs of branch cracks are not always proportional to √r. The energy release rate is obtained from the measured CODs, and it is found that energy release rate is continuous at bifurcation point in both cases of Homalite 100 and Araldite B. Energy flux that shows the energy flow toward a crack tip is also obtained.     

A DECAGONAL QUASICRYSTAL WITH AN ARC-SHAPED CRACK

A decagonal quasicrystal,which is weakened by an arc-shaped crack penetratingthrough the solid in the period direction,and which is subjected to remote uniform phonon stresses,is investigated by applying the complex variable method which is just developed by the authors.It is found that the phonon and phason stresses near the crack tips exhibit inverse square rootsingularities.The four complex stress functions characterizing the phonon and phason fields arederived.Explicit expressions for the phonon and phason stress intensity factors,crack openingdisplacements and energy release rate are also presented.     

Energy release rates in a constrained epoxy disc with Hookean and Mooney–Rivlin materials

Assuming that the disc material can be modeled either as Mooney–Rivlin or as Hookean and the steel ring enclosing the disc as Hookean, the energy release rates as a function of the crack length are evaluated and compared. Two loadings are considered––one in which the surface of the star shape hole in the disc is loaded by a uniform pressure and the other in which the temperature of the composite body is uniformly raised. It is found that the linear and the nonlinear analyses give qualitatively similar results for the two loadings. For each load, the energy release rate increases with an increase in the starter crack length, reaches a maximum value and then decreases gradually.     

Shear Failure of Inconel 718 under Dynamic Loads

Kinetics of deformation and fracture of nickel–iron alloy Inconel 718 under dynamic shear loading was measured using a split torsional Hopkinson bar facility and high-speed photography. Tubular specimens with a reduced gage length and a starter notch were sheared at strain rates up to 6 × 10³ s⁻¹. High-speed photographs of fiducial lines scribed on the specimen surface showed the development of local strains and cracking. This paper describes the experimental and analytical procedures, illustrates average and local plastic strain evolution, and presents shear crack initiation times and propagation speeds.     

Drying of porous oil shales

Rate processes including change of phase are modelled analytically for a half-space porous substance exposed to a jump in external temperature. The model predicts, in a closed form, the pressure build-up and the rate of evaporation of volatiles from a porous matrix. It assumes two distinct regions separated by a moving interface where the change of phase takes place. One region maintains its initial concentration of volatiles while the second is devoid of volatiles. Different thermophysical properties are considered for the two regions. The model was applied for the evaporation of moisture from oil-shale. Results are given in a parametric form     

船用气水分离器惯性级流场分析及阻力特性研究

采用二阶全展开ETG有限元与大涡模拟（LES）相结合的算法，对120－20－35－3型船用气水分离器惯性级在不同雷诺数条下的流动进行了模拟，通过其中一组雷诺数条件下计算所得该实验件阻力系数与物理实验所得阻力系数相比较，确定出该雷诺数条件下采用大涡模拟时所需的亚格子应力模型常数，将该常数带入其它各组雷诺数条件下的计算中，并将计算结果与相同条件下的物理实验结果相比较，证实了该常数的通用性。该常数一经确定，对各雷诺数条件下的流场进行分析，结果反映出采用二阶全展开ETG有限元与大涡模拟（LES）相结合的算法可以捕捉到非常丰富的涡系及涡动的时变过程。在流场分析的基础上本文计算了该实验件内的能耗场，计算结果表明实验件内的能耗主要集中在大涡丰富的区段内。